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Abstract. The constraints obtained by several authors (Eke et al. 1998; Henry 2000) 

^ ' on the estimated values of f2 m , n and as are revisited in the light of recent theoretical 
ID [ 

qq I developments: 1) new theoretical mass functions; 2) a more accurate mass-temperature 
relation, also determined for arbitrary f2 m and Qa- I re-derive the X-ray Temperature 
Function (XTF), similarly to Henry (2000) and Eke et al. (1999), re-obtaining the con- 
' straints on Q m , n, as- The result shows that changes in the mass function and M-T 
' relation produces an increase in Q m of ~ 20% and similar results in erg and n. 
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1. Introduction 

Galaxy clusters represents the virialization stage of exceptionally high peaks of initial 
density perturbations on comoving scales of ~ 10/i -1 Mpc, and as such they provide useful 
cosmological probes. The evolution in the abundance of clusters is strongly dependent 
on the cosmological density parameter, Q m (Evrard 1989; Oukbir & Blanchard 1992; 
Eke et al. 1996; Donahue et al. 1998; Borgani et al. 1999). X-ray observations provide 
a very efficient method to identify distant clusters down to a given X-ray flux limit, 
and hence within a known survey volume for each luminosity, Lx- For this reason, most 
studies using clusters as cosmological probes are based on X-ray samples. It is well known 
that clusters are strong X-ray emitters and so cluster evolution can be inferred from the 
study of X-ray properties of distant clusters. The amount of observational data concerning 
high-redshift cluster properties has increased in the past years. EMSS (Henry et al. 1992; 
Gioia & Luppino 1994), ASCA (Donahue 1996; Henry 1997) and ROSAT (Ebeling et 
al. 1997; Rosati et al. 1998) In addition, galaxy velocity dispersions for a well-defined 
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sample of high-redshift clusters (Carlberg et al. 1996) are provided by the CNOC survey. 
The Press & Schechter (1974) (hereafter PS) formalism has been heavily used to model 
the cluster population. The combination of the PS mass function and the X-ray cluster 
catalogs represents a unique opportunity to constraint cosmological parameters, (e.g. 
the mass density parameter, Q m ). Although the analytical framework of the PS model 
has been greatly refined and extended (e.g., Lacey & Cole 1993), it is well known that 
the PS mass function, while qualitatively correct, disagrees with the results of N-body 
simulations. The quoted discrepancy is not surprising since the PS model, as any other 
analytical model, should make several assumptions to get simple analytical predictions. 
There are different methods to trace the evolution of the cluster number density: 

a) The X-ray temperature function (XTF) has been presented for local (e.g. Henry & 
Arnaud 1991) and distant clusters (Eke et al. 1998; Henry 2000). The mild evolution of 
the XTF has been interpreted as a strong indication for a low density universe (0.2 < 
O m < 0.6). 

b) The evolution of the X-ray luminosity function (XLF). 

The results for f2 m obtained span the entire range of acceptable solutions: 0.2 < Q m < 
1 (see Reichart et al. 1999). The reasons leading to the quoted discrepancies has been 
studied in several papers (Eke et al. 1998; Reichart et al. 1999; Donahue & Voit 1999; 
Borgani et al. 2001). 

Although the quoted uncertainties has been so far of minor importance with respect 
to the paucity of observational data, a breakthrough is needed in the quality of the 
theoretical framework if high-redshift clusters are to take part in the high-precision-era 
of observational cosmology. 

These reasons lead me to re-calculate the constraints on il m , n and er 8 , using the 
XTF. In Sect. 2, I re-calculate the XTF, as done by Henry (2000) and Eke et al. (1998) 
and obtained constraints for Q m and n and cr 8 . Sect. 3 is devoted to results and to 
conclusions. 

2. Constraints to cosmological parameters from the XTF 

The mass function (MF) is a critical ingredient in putting strong constraints on cosmolog- 
ical parameters (e.g., fi m ). Observationally the local mass function has been derived from 
measuring masses of individual clusters from galaxy velocity dispersions or other optical 
properties by Bahcall and Cen (1993), Biviano et al. (1993), and Girardi et al. (1998). 
However, the estimated virial masses for individual clusters depend rather strongly on 
model assumptions. As argued by Evrard et al. (1996) on the basis of hydrodynamical 
N-body simulations, cluster masses may be presently more accurately determined from 
a temperature measurement and a mass-temperature relation determined from detailed 




Fig. 1. (a) Comparison of various mass functions. The dotted line represents Sheth & 
Tormen (2002) prediction, the solid line that of Jenkins et al. (2001) and the dashed line 
that of Del Popolo (2000b). 



observations or numerical modeling. Thus alternatively, as a well-defined observational 
quantity, the X-ray temperature function (XTF) has been measured, which can be con- 
verted to the MF by means of the mass-temperature relation. 

The cluster temperature function is defined as: 

N(T,z) = N(M,z) — (1) 

While the mass function, N(M,z), gives the mass and distribution of a population of 
evolving clusters, the Jacobian describes the physical properties of the single cluster. 

Comparison of the predictions of the PS theory with the SCDM and OCDM cosmolo- 
gies, performed by Tozzi & Governato (1998) and Governato et al. (1999), have shown 
discrepancies between PS predictions and N-body simulations, increasing with increasing 
z. 
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(c) (d) 

Fig. 2. (a) The 68% confidence contours for the parameters n, k Q and Q m for the open 
model. The dashed lines arc lines of constant a%. (b) A (likelihood) for the parameter 
fi m . The solid line is obtained from the model of this paper while the dotted line is that 
calculated by Henry (2000). The dashed lines represent various confidence levels (65%, 
90%, 95%, 99%). (c)The 68% confidence contours for the parameters ag, and f2 m for the 
open model (see also Henry (2000), Fig. 9). (d) Constraints on J7a and f2 m obtained using 
the same 25 clusters used in Henry (2000), for the local sample, while the high rcdshift 
sample is constituted from all the EMSS clusters with z > 0.3 and RX J0152. 7-1357 (see 
Henry 2002). The solid lines are the 1 and 2 a contours obtained using the mass function 
and the M-T relation of this paper, while the dashed line is the 1 a contour obtained 
using the PS mass function and the M-T relation of Pierpaoli et al. 2001. 



In the following, I'll use the mass function modified as described in Del Popolo (2000b) 
and an improved form of the M-T relation in order to calculate the mass function (Del 
Popolo 2000a). These are respectively given by: 



n{m, z) ~ 1.21 
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where a = 0.707. and 
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(see Del Popolo 2000a for a derivation), where M ~ 5 x 10 14 /i 1 M Q , i n = 

mi = 5/(n + 3) (being n the spectral index), and: 

H (i-a m -n A )S 

K(mi,x) — Fx (mi — 1) LerchPhi(x, 1, 3mi/5 + 1) — 



where f 1 is defined in Del Popolo (2000a, b) and the LerchPhi function is defined as 
follows: 



where x = 1 + (^-) 2 / 3 which is connected to mass by M = M x~ 3 ™ 1 / 5 (V2000), and 
where Ko(mi,x) indicates that K(mi,x) must be calculated assuming t = tQ. Since Eq. 
J2J and Eq. © where obtained in previous papers, Del Popolo 2002a, b, and since there 
is a wide description of the underlying assumptions, I refer the reader to those papers 
for details. 

Before going on, I want to discuss about the use of Eq. (|2J). Jenkins et al. (2001), 
obtained a mass function, which is regarded as perhaps the most accurate model to date, 
from the Hubble volume simulations. Jenkins et al. (2001) showed that, although the 
mass functions in their simulations scaled in accordance with the excursion set prediction, 
Sheth & Tormen 2002 (Eq. (2)) slightly overestimated the unconditional mass functions 
in their simulations. It interesting to note that as shown by Pierpaoli et al. (2001), 
using Sheth & Tormen (2002) or Jenkins et al. (2001, 2003) mass function, rather than 
PS, the value of parameters like eg are changed by a small amount (4-8 %). Sheth & 
Tormen (2002) also showed that changing the parameter a in their Eq. (2) from 0.707 
to 0.75 reduces the discrepancy between it and the simulations substantially. In Del 
Popolo (2002a), I showed that Eq. of this paper, in agreement with Jenkins et al. 
(2001), predicts smaller values of the mass function expecially for high v, with respect 
with Sheth & Tormen's predictions. In other terms, Eq. ((2J of this paper is in very 
good agreement with Jenkins et al. (2001). In Fig. 1 I plot a comparison of the various 
mass functions: the dotted line represents Sheth & Tormen (2002) prediction, the solid 
line that of Jenkins et al. (2001) and the dashed line that of Del Popolo (2000b). As 
shown, Jenkins et al. (2001) mass function is almost indistinguishable from that used 
in this paper. For what concerns small masses, our formula, in agreement with Sheth & 
Tormen 2002, differs dramatically from the one Jenkins et al. (2001), propose. Simulations 
currently available do not probe the regime where v < 0.3 or so (the Jeans mass is at 



F (mi — 1) LerchPhi(x, 1, 3mi/5) 



(4) 




(5) 
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about v < 0.03 (Sheth & Tormen 2002)). New simulations are needed to address which 
low mass behavior is correct. In other terms, the mass function obtained in this paper is 
in very good agreement with Jenkins et al. (2001) in the regime probed by simulations. 
Moreover, the constraints obtained in almost all the papers in literature used the PS mass 
function except a few papers (e.g. Borgani et al. 2001; Henry 2002) and the M-T relation 
the usual one obtained from the virial theorem. In other words, this paper introduces 
noteworthy improvements on the previous calculations in literature. The mass entering in 
Eq. (J2J) should be interpreted as the mass contained inside a radius, riso encompassing a 
mean overdensity p = 180/5. However, scaling relations connecting mass X-ray observable 
quantities may provide the mass at different values of p/p. In this case we follow White 
(2001) and rescale the masses assuming an NFW (Navarro, Frenk & White 1996) profile 
for the dark matter halo with a concentration c = 5 appropriate for rich clusters (see 
also Pierpaoli et al. 2001, 2003; Schuecker et al. 2002 for more details). 

Similarly, it necessary to give arguments to use Eq. instead of the new mass/X-ray 
temperature relations obtained from simulations or Chandra data within the last year 
(see, e.g., Pierpaoli et al. 2001, 2003 for a reference). As shown in Del Popolo (2002a), 
Eq. J2J reduces to a similar equation to that used in Pierpaoli et al. (2001) (Eq. 13), in 
the early-time, namely: 



where E(z) 2 — f2jyr(l + z ) 3 + ^A, and the term depending on fl\ in Eq. (13) of Pierpaoli 
et al. (2001) (which is a correction to the virial relation arising from the additional r 2 
potential in the presence of A) is neglected since it produces only a small correction (see 
Pierpaoli et al. 2001). I also want to add that Eq. (13) of Pierpaoli et al. (2001) or Eq. (4) 
of Pierpaoli et al. (2003), comes from rather simplistic arguments (dimensional analysis 
and an assumption that clusters are self-similar) and is a good approximation to both 
observations and simulations but this last are sufficiently computationally demanding 
that they cannot explore parameter space efficiently and so it is necessary to determine 
coefficients by means of simulations, while scalings are taken from simple theoretical 
models (Pierpaoli et al. 2001). Eq. (13) of Pierpaoli et al. (2001), is valid for systems 
hotter than about 3 keV. We know that recent studies have shown that the self-similarity 
in the M-T relation seems to break at some keV (Nevalanien et al. 2000; Xu, Jin & Wu 
2001): Afshordi & Cen (2001) has shown that non-sphericity introduces an asymmetric, 
mass dependent, scatter for the M-T relation altering its slope at the low mass end (T ~ 3 
keV). These effects are taken into account by mine M-T relation, which as previously 
told gives same results of that of Pierpaoli et al. (2001, 2003) in the range of energy in 
which this is valid. For what concerns the definition of masses in the M-T relation, the 
mass is defined, according to Del Popolo (2000a), Voit & Donahue 1998 and Voit 2000, 



M oc T 3 /V 1/2 A C ~ 1/2 oc T^ 2 (A C E 2 )-^ 2 



(6) 
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(a) (b) 



Fig. 3. A(— lnlikclihood) for f2 m (left panel) and as (right panel), marginalized over the 
other two parameters ((r, erg), left panel and (O m , T), right panel). Solid lines are the 
prediction obtained from the model of this paper while the dotted lines those obtained 
from Eke et al. (1998). Dashed lines are 1 and 2 a significance and the 3 a corresponds 
to the top of each panel. 



as M2oo- In order to combine the mass function and M-T relation, I performed a mass 
transformation, as above reported (see also Pierpaoli et al. 2001, 2003; Schuecker et al. 
2002 for more details). 

Introducing Eq. (J2J into Eq. JJ) and using the M-T relation, in the peculiar case that 
the variance is given by: 



a = C( ,f )-( 3 +">/ 6 = Cm-^' & 
y 10 15 M & ' 

where n is the spectral index. 
I get: 

dN _ N\ (m, z) 
~dT ~ N 2 (m,z) 



(7) 



(8) 
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K(rn± ,x) 



— + t? + 



(j\/ 15 m)S/3 



_L , j.3 I K (rni,x) 



(9) 



where, I have defined: m = j^, mi = 5/(n + 3), m 2 = "-/6, m 3 = 1 + ra/3, rn 4 = 
-l/2-n/6, a, a x = 0.06, and a 2 = 0.57 a 3 = -5/(3ml), my = 3ml/5 + l, £i = (io/i) 2/3 , 
r = 2/3, g = 8/3, c5 c0 (z) = 5, and LerchPhi(m a3 , 1,toi;) = LP/i. 

In order to use the same notation and variance of Henry (2000), the constant C is 
defined as: 

O = 0.675 / r(3 + n) 8 in (1 /2 r )- / 

y 2™n(2 + n)(l-n)(3-n) V / 

and then erg = cr(il m , M = 0.594 x 10 15 /i _1 il m ), and the XTF depends on the parameters 
n, fc Q and Jl m . The data that shall be fitted to the theory previously described, are 
those described in Section. 2 of Henry (2000). I use a maximum likelihood fit to the 
unbinned data in order to determine various model parameters. The method is described 
in Marshall et al (1983). Using Poisson probabilities, we have: 

N 

Li = Y[[X(z t ,Li)dzdLe- x{z ^ )dzdL } J[ e ~K^,L 3 )dzdL ^ 

i 3 

where \(z, L)dzdL = p(z, L)Q(z, L)^dzdL is the expected number of objects in dzdL 
at z, L in the sample and the index j runs over all the differential elements in which no 
objects were observed. Defining S = —2ln{Li) and dropping terms independent of the 
model parameters one finds: 

^ = -2^1np(z 2 ,L. i ) + 2 f f p( z ,L)n(z,L) ^^ dzdL (12) 

Following (Henry 2000) notation and assumption Eq. I|12l) is written as: 

;2i 



2Y> La^z^kT ^^^ } +2 r^dkT 
~[ l dzdn J A-T min 



'" n(Q m , z, kT)Q(z, kT) Z) dz (13) 

(Henry 2000), where N is the number of clusters observed, n(f2 m , z, kT) is the temperature 
function, Q(z, kT) is the solid angle in which a cluster with temperature kT at redshift z 
could have been detected (the selection function) and d V j^ l ' Zl ' > is the differential volume, 
which is given in the Appendix of Henry (2000) (see also Henry 2000 for a description 
of d V }2m' Z ^ m t nc fi rs t term of Eq. I|13|l). The best estimates for the parameter are 
obtained minimizing S. 

At this point, we can fit the data described in Section. 2 of Henry (2000) to the theory 
previously described using the quoted maximum likelihood method. The most general 
description of the results requires the three parameters of the fit. I show these results in 
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Fig. 2, where I plotted the results for the open model. It is straightforward to read off the 
value of n, which is at the 68% confidence for the open model. These values shows that 
the correction introduced by the new form of the mass function and M-T relation gives 
rise to higher values of ft m (Sl m = 0.6±0.13, while it is f7 m = 0.49±0.12 for Henry (2000)) 
and n = -1.5 ± 0.32 (n = -1.72 ± 0.34 in Henry (2000)). The presentation in Fig. 2a 
is somewhat difficult to appreciate, so we also give the constraints for fewer parameters. 
In Fig. 2b, I plot A(likelihood) for the parameter f2 m . The solid line is obtained from 
the model of this paper while the dotted line is that calculated by Henry (2000). The 
dashed lines represent various confidence levels (68%, 90%, 95%, 99%). Constraints are 
relatively tight when considering this single parameter. I find that f2 m = 0.6^o 11 a t the 
68% confidence level and at £l m = 0.6^q'2 3 the at 95% confidence level for the open 
model. 

The constraints in Fig. 2a are plotted into a more conventional format in Fig. 2c. 
Three parameters are still required, but the constraints on n and kg are collapsed into 
og- Fig. 2c, plots the 68% confidence contours for the parameters a$, and Sl m for the 
open model (see also Henry (2000), Fig. 9). In Fig. 2d, I plot the constraints on 12 a and 
fi m obtained using the same 25 clusters used in Henry (2000), for the local sample, while 
the high redshift sample is constituted from all the EMSS clusters with z > 0.3 and RX 
J0152. 7-1357 (see Henry 2002). The solid lines are the 1 and 2 a contours obtained using 
the mass function and the M-T relation of this paper, while the dashed line is the 1 a 
contour obtained using the PS mass function and the M-T relation of Pierpaoli et al. 
2001. 

For a CDM spectrum, the expression for the XTF is much more complicated. It can 
be obtained combining Eq. J2J) , Eq. Q , and our M-T relation. The mass variance can be 
obtained once a spectrum, P(k), is fixed, by: 



1 f°° 

a 2 (M) = \ dkk 2 P(k)W' 
2?H J 



(kR) (14) 



where W(kR) is a top-hat smoothing function: 

W{kR) = -^-^ (sinfci?- kRcoskR) (15) 

and the power spectrum P(k) = Ak n T 2 (k) is fixed giving the transfer function T(k). 
The CDM spectrum used in this paper is that of Bardeen et al. (1986) (equation (G3)). 
The shape of the spectrum and its amplitude at 8 hr x Mpc are described by as and T, 
respectively. Similarly to Eke et al. (1998) the fitting parameters were f2 m , T, as, but 
just the results for fl m and erg were plotted. Using the data used by Eke et al. (1998) 
(see their Sec. 3 1 ) and the maximum likelihood parameter estimation given in Eke et al. 

1 They combined the temperature data for 25 local clusters by Henry & Arnaud (1991) with 
the sample of 10 EMSS clusters at 0.3 < z < 0.4 by Henry et al. (1997) 
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(1998), Sec. 4, I obtain the results plotted in Fig. 3a-b. Fig. 3a plots the A(-lnlikelihood) 
for fi m (left panel) and erg (right panel), marginalized over the other two parameters ((r, 
as) left panel, and (fi m , T), right panel). Solid lines are the prediction obtained from the 
model of this paper while the dotted lines those obtained from Eke et al. (1998). Dashed 
lines are 1 and 2 a significance and the 3 a corresponds to the top of each panel. Fig. 
3 shows that £! m (left panel) and a% (right panel) are increased with respect to Eke et 
al. (1998) prediction: while in Eke et al. (1998) fl m = 0.52±°;^, and a 8 = 0.63±g;g|, I 
find that fi m = 0.62+o'*5 and as = 0.731q;qq. This shows again an increase in O m , also 
in qualitative agreement with Eke et al. (1998) calculation taking account of changes in 
the threshold for collapse suggested by Tozzi & Governato (1998). 



3. Results and discussion 

In this paper, I have revisited the constraints obtained by several authors (Eke et al. 
1998; Henry 2000) on the estimated values of O m , n and trg in the light of recent theo- 
retical developments: new theoretical mass functions, a more accurate mass-temperature 
relation, also determined for arbitrary fl m and Qa- I repeated the Henry (2000) analysis 
but differently from the quoted paper, I changed the mass function and M-T relation, 
adopting again those of Del Popolo (2000a, b). The new form of the mass function and 
M-T relation gives rise to higher values of fi m (f2 m = 0.6 ± 0.13 in my estimation, while 
it is fl m = 0.49 ± 0.12 for Henry (2000)) and n = -1.5 ± 0.32, in my estimation, while 
n = —1.72 ± 0.34 in Henry (2000). I have also obtained some constraints on Qa and f2 m 
obtained using the same 25 clusters used in Henry (2000), for the local sample, while 
the high redshift sample is constituted from all the EMSS clusters with z > 0.3 and 
RX J0152. 7-1357 (see Henry 2002). The 1 a contours obtained using the mass function 
and the M-T relation of this paper shows that for A = 0, it is 0.32 < £! m < 0.57 in 
the case of Henry (2002) and 0.43 < £! m < 0.73 in my estimation. The figure shows the 
constraints to f2 m for different values of Similar results to that obtained in the pre- 
vious comparison with the Henry (2000) results are obtained changing data and method, 
by following Eke et al. (1998). I obtain a value of £l m = 0.621^;^, while in Eke et al. 
(1998) ft m = 0.52+° 6 17 , and a 8 = 0.731^6 whilc in Ekc et aL ( 1998 ) °8 = 0.63±g;g|. 
This shows again an increase in fi m , also in agreement with Eke et al. (1998) calculation 
taking account of changes in the threshold for collapse suggested by Tozzi & Governato 
(1998). 2 As previously told, this paper has its aim that of studying how "systematic 
uncertainties" can influence the values of some cosmological parameters. It is well known 



2 Replacing S c with the <5 c ff of Tozzi & Governato (1998) , produces qualitatively similar changes 
on the mass function as those obtained using the Del Popolo (2002b) mass function. 
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that in literature the values obtained for f2 m span the range 0.2 < O m < 1 (Reichart et 
al. 1999). 

Sadat et al. (1998) and Reichart et al. (1999) Blanchard & Bartlett (1998) found 
results consistent with f2 m = 1. Viana & Liddle (1999) found that fi m = 0.75 with 
O m > 0.3 at the 90% confidence level and Q m ~ 1 still viable. Blanchard et al. (1998) 
found almost identical results (f2 m ~ 0.74 with 0.3 < £! m < 1.2 at the 95% confidence 
level). Eke et al. (1998) found Q m = 0.45 ± 0.2. It is interesting to note (as previously 
mentioned) that Viana & Liddle (1999) used the same data set as Eke et al. (1998) and 
showed that uncertainties both in fitting local data and in the theoretical modeling could 
significantly change the final results: they found fl m ~ 0.75 as a preferred value with a 
critical density model acceptable at < 90% c.l. 

Different results were obtained by Bahcall et al. (1997) (fl m — 0.3 ± 0.1), Fan et al. 
(1997) (O m = 0.3±0.1), Bahcall & Fan (1998) (O m = 0.2±g f ) and several other authors. 
More recent studies has been published, giving constraints on cr 8 and fl m : Pierpaoli et 
al. (2003) using a sample updated from HIFLUGCS they found a 68 % confidence range 
of o- 8 = 0.77t° ° o l for a standard ACDM model, while using the XLF from the REFLEX 
survey they obtain cr 8 = 0. 79^^07 f° r tnc same standard ACDM model. Schuecker et 
al. (2002), using the REFLEX cluster sample and analyzing simultaneously the large 
scale clustering and the mean abundance of galaxy clusters get precise constraints on 
the normalized cosmic matter density fi m (fi m = 0.34lto;o2g) an d the linear theory 
r.m.s. fluctuations in mass, as (&s = 0. 711+0. 039_o.03i)- WMAP (Wilkinson Microwave 
Anisotropy Probing) precision data enabled accurate testing of cosmological models. 
Spergel et al. (2003) found that Qh 2 = 0.135t o ;{$g, a s = 0.84±0.04, where h = 0.71±g;gf . 
The formal random errors given in WMAP are of the 5 % level while they are higher in our 
study. At this regard, I want to recall that the principal goal of the present study is not 
that of determining with high precision the values of the quoted cosmological parameters 
but to show how different values 3 , discrepancies in the values of cosmological parameter 
quoted in literature are connected to a different mass function and M-T relation (the 
basic aim of the present work is to study how systematic uncertainties in the application 
of the X-ray temperature function and mass/X-ray temperature relation can influence 
the estimates of important cosmological parameters). 

The reasons leading to the quoted discrepancies has been studied in several papers 
(Eke et al. 1998; Reichart et al. 1999; Donahue & Voit 1999; Borgani et al. 2001). Ac- 
cording to Reichart (1999) unknown systematic effects may be plaguing great part of 
the quoted results. Our analysis shows that improvements in the mass function and M-T 

3 As told in the introduction the results for f2 m obtained span the entire range of acceptable 
solutions: 0.2 < f2 m < 1 
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relation increases the value of £l m . The effect of this increase is unable to enhance signif- 
icantly the probability that f2 ro = 1 in the case of constraints like that of Fan, Bahcall & 
Cen (1997) (fi m = 0.3) or Bahcall & Fan (1998) (Q m = 0.2), and can give a small contri- 
bution even in the case of larger values for the value of the constraints obtained. However, 
in any case it shows that even small correction in the physics of the collapse can induce 
noteworthy effects on the constraints obtained. Moreover, even changing the data or the 
way they are analyzed gives different results. As an example, changing their low-redshift 
sample, Donahue & Voit (1999) showed a change in Q m from 0.45 to 0.3. Furthermore, as 
observations are reaching the first epoch of cluster assembly, treating them as dynamical 
relaxed and virialized systems is undoubtly an oversemplification. Hierarchical clustering 
scenario predicts that a fraction between 0.3 and 0.6 of the z = 1 population of clusters 
are observed less than 1 Gyr after the last major merger event and then are likely to be 
in a state of non-equilibrium. Although the quoted uncertainties has been so far of minor 
importance with respect to the paucity of observational data, a breakthrough is needed 
in the quality of the theoretical framework if high-redshift clusters are to take part in 
the high-precision-era of observational cosmology. 
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